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The Ra and Pr number scaling of the Nusselt number Nu, the Reynolds number Re, the tem- 
perature fluctuations, and the kinetic and thermal dissipation rates is studied for (numerical) ho- 
mogeneous Rayleigh-Benard turbulence, i.e., Rayleigh-Benard turbulence with periodic boundary 
conditions in all directions and a volume forcing of the temperature field by a mean gradient. This 
system serves as model system for the bulk of Rayleigh-Benard flow and therefore as model for 
the so called "ultimate regime of thermal convection". With respect to the Ra dependence of Nu 
and Re we confirm our earlier results Q which are consistent with the Kraichnan theory and 
the Grossmann-Lohse (GL) theory HUSH, which both predict Nu ~ Ra 1/2 and Re ~ Ra 1/2 . 
However the Pr dependence within these two theories is different. Here we show that the numerical 
data are consistent with the GL theory Nu ~ Pr 1 / 2 , Re ~ Pr -1 ^ 2 . For the thermal and kinetic 
dissipation rates we find ee/(K,A 2 L- 2 ) ~ (RePr)°- S7 and e u /(!/ 3 L~ 4 ) ~ Re 2 77 , also both consistent 
with the GL theory, whereas the temperature fluctuations do not depend on Ra and Pr. Finally, 
the dynamics of the heat transport is studied and put into the context of a recent theoretical finding 
by Doering et al. ■ 
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The scaling of large Rayleigh number (Ra) Rayleigh- 
Benard convection has attracted tremendous attention in 
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the last two decades 

3 El El El E] 

H Ull Elli lHal^riHllHJ There is increasing 
agreement that in general there are no clean scaling laws 
for Nu(Ra, Pr) and Re(Ra, Pr), apart from asymptotic 
cases. One of these asymptotic cases has been doped the 
"ultimate state of thermal convection" , where the heat 
flux becomes independent of the kinetic viscosity v and 
the thermal diffusivity k. The physics of this regime is 
that the thermal and kinetic boundary layers have broken 
down or do not play a role any more for the heat flux and 
the flow is bulk dominated. The original scaling laws 
suggested for this regime are 



Nu ~ Ra 1/2 (log Ra)- 3/2 Pr 1/2 
Re ~ Ra 1/2 (log Ra)- 1/2 Pr~ 1/2 , 
for Pr < 0.15, while for 0.15 < Pr < 1: 

Nu ~ Ra 1/2 (log Ra)- 3/2 Pr- 1/4 
Re ~ Ra 1/2 (log Ra)-V 2 Pr- 3 ? 4 . 



(1) 
(2) 

(3) 
(4) 



The GL theory also gives such an asymptotic regime 
which is bulk dominated and where the plumes do not 
play a role (regimes IV/ and IVJ of Refs. H&HH). 
It has the same Ra dependence as in eqs. but 
different Pr dependence, namely 



Nu ~ Ra 1/2 Pr 1/2 , 



As a model of the ultimate regime we had suggested 
.1] homogeneous RB turbulence, i.e., RB turbulence with 
periodic boundary conditions in all directions and a vol- 
ume forcing of the temperature field by a mean gradient 

M, 



do 



(u-8)0 = kO 2 



A 
— 'i 
L 



(7) 



Here = T + (A/L)z is the deviation of the tem- 
perature from the linear temperature profile —(A/L)z. 
The velocity field u(x, t) obeys the standard Boussinesq 
equation, 



c)u 
~dt 



(it ■ d)u = - Vp + ud 2 u + (3gz9. 



(8) 



Here, f3 is the thermal expansion coefficient, g gravity, 
p the pressure, and 0(x,t) and Ui(x,t) are temperature 
and velocity field, respectively. Indeed, in Ref. [| we 
showed that the numerical results from eqs. @ and JSJ) 
are consistent with the suggested 0, 0, S Wm Ra de- 
pendence of Nu and Re, Nu ~ Ra 1 / 2 and Re ~ Ra 1 / 2 . 
However, the Pr dependences of Nu and i?e, for which 
the predictions of Kraichnan and GL 0, S 0. @ are 
different, has not yet been tested for homogeneous tur- 
bulence: this is the first aim of this paper (Section III). 
Section II contains details of the numerics. In Section IV 
we study the bulk scaling laws for the thermal and kinetic 
dissipation rates and compare them with the GL theory. 
In that Section we study the temperature fluctuations 

1/2 

(5) 9 1 = (9 2 ) '-. The dy namics of the flow, including Nu(t) 
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and its PDF (probability density function), is studied in 
Section V and put into the context of a recent analyt- 
ical finding by Doering and coworkers Q. Section VI 
contains our conclusions. 



II. DETAILS OF THE NUMERICS 

Our numerical simulation is based on a Lattice Boltz- 
man Equation (LBE) algorithm on a cubic 240 3 grid. 
The same scheme and resolution has already been used 
m We run two sets of simulations in statistically 

stationary conditions. The first at fixed Pr = 1 varying 
the Ra number between 9.6 ■ 10 4 and 1.4 • 10 7 . The sec- 
ond at fixed Ra = 1.4 • 10 7 . This, the highest value we 
can reach at the present resolution, was studied for five 
different Pr numbers, 1/10, 1/3, 1, 3 and 4. We recorded 
shortly-spaced time series of Nu and root mean squared 
(rms) values of temperature and velocity and we stored a 
collection of the whole field configurations, with a coarse 
time-spacing. The length of each different run ranges 
between 64 and 166 eddy turn over times. Our simu- 
lation was performed on a APEmil/e machine in a 128 
processor configuration [5£|- Each eddy turnover 

times requires on average 4 hours of computation. The 
total computational time required for the whole set sim- 
ulations is roughly 150 days. The total number of stored 
configurations is around 2000. 



III. Nu{Ra,Pr) AND Re{Ra,Pr) 

The Nusselt number is defined as the dimensionless 
heat flux 

Nu = (<u 3 T) A>t (z) - k (d 3 T) A t (zj) 

(u 3 0) A t (z) 



kAL- 
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(9) 



where the average {■■■) At is over a horizontal plane and 
over time. From eqs. (0-© one can derive two exact re- 
lations for the volume averaged thermal dissipation rate 
eg = n((di0) 2 ) v and the volume averaged kinetic dissi- 
pation rate e u = v ((SiWj) 2 ) „, namely 



e u = — rNuRaPr 
L 4 



A 2 
k-—Nu. 



(10) 



(11) 



One can therefore numerically compute Nu in three 
different ways: (i) from its direct definition 0, (ii) from 
the volume averaged kinetic dissipation rate (|10|) . (iii) 
from the volume averaged thermal dissipation rate 111(1 . 

The results are shown in Figure^, as a function of Ra 
for Pr = 1. There is very good agreement of Nu obtained 



from the three different methods for all Ra, giving us 
further confidence in the convergence of the numerics. If 
we fit all data points beyond Ra = 10 5 with an effective 
power law, we obtain Ra ~ i?a°- 50±0 05 , consistent with 
the asymptotically expected law Nu ~ Ra 1 / 2 [53 

In Figure ^ we display Nu as function of Pr for fixed 
Ra = 1.4 • 10 7 . For the cases with Pr ^ 1 the conver- 
gence of the three different methods to calculate Nu is 
not perfect. This may be due to numerical errors in the 
resolution of the small scale differences, especially when 
v and k are considerably different. However, one can 
clearly notice a strong increase of Nu with Pr. A fit with 
an effective power law gives Nu ~ p r °- 43 ± 07 ; which is 
consistent with the asymptotic power law Nu ~ Pr 1 / 2 
suggested by the GL theory and by the small Pr regime 
suggested by Kraichnan, but not with Kraichnan's 
large Pr regime (J3J). Increasing further Pr (at fixed Ra) 
the flow will eventually laminarize, i.e., can no longer be 
considered as model system for the bulk of turbulence. 
This also follows from Figure \Bp 1 in which we show the 
Reynolds number: 



Re = 



i'L 



(12) 



as function of Pr for fixed Ra — 1.4- 10 7 . Note that this 

is the fluctuation Reynolds number, defined by the rms 

1/2 

velocity fluctuation u' = (it 2 ) : in homogeneous RB 
no large scale wind exists. Re(Pr) displays an effective 
scaling law Re ~ p r -o.55±o.oi^ cons i s t e nt with the GL 
prediction Pr~ 1//2 for the ultimate regime (if one identi- 
fies the wind Reynolds number in GL with the fluctuation 
Reynolds number here) and also with the Kraichnan pre- 
diction (JIJ . Also the Ra scaling of Re is consistent with 
GL (and also with Kraichnan), Re ~ Ra 1 / 2 , as seen from 
Figure and as already shown in Ref. 0. 



IV. SCALING LAWS FOR e„, e AND THE 
TEMPERATURE FLUCTUATIONS 



A. Kinetic and thermal dissipations 

The homogeneous RB turbulence offers the opportu- 
nity to numerically test one of the basic assumptions of 
the GL theory, namely, that the energy dissipation rate 
in the bulk scales like 



L 4 



Re 6 



(13) 



In Figure |3fa) we plot e u /(v 3 L 4 ) vs. Re for all Ra 
and Pr and find e u /(v 3 L.- 4 ) - i?e 2 - 77±0 03 , close to the 
expectation itl^l . 

The disentanglement of the thermal dissipation rate eg 
into two different scaling contributions is less straightfor- 
ward. The GL theory decombincs: 



eg = c 3 {RePr) 1/2 + c 4 (RePr) 



(14) 
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FIG. 1: (a) Nu(Ra) for Pr = 1, computed in three different 
ways: (•) using Eqn. ©, (□) using Eqn. (IIOL and (o) from 
Eqn. (jnj- The power law fits, performed on the mean value 
of the three different estimates and for Ra > 10 s , gives a slope 
0.50 ± 0.05. (b) Nu(Pr) for Ra = 1.4 • 10 7 , fit performed as 
before, with a resulting slope of 0.43 ± 0.07. 



FIG. 2: (a) Re(Ra) for Pr = 1, with a fitted slope 0.50±0.02. 
(b) Re(Pr) for Ra = 1.4- 10 7 , with a fitted slope -0.55±0.01. 



B. Temperature fluctuations 



where the first term has been interpreted as boundary 
layer and plume contribution eg w and the second one 
as background contribution ee,6g p|- The prefactors C3 
and C4 are given in Ref.|jL Plumes are interpreted as 
detached boundary layer Q. For homogeneous RB tur- 
bulence one would expect the background contributions 
to be dominant as there is no boundary layer. But 
still some plumes may also develop in the bulk and this 
is confirmed by the fact that we find a scaling law in 
between the asymptotes {RePr) 1 / 2 and RePr, namely, 
ee ~ (RePr) ' 87 04 : closer to the background behavior 
just as one would guessed. 



In our numerics we find the temperature fluctuations 
1/2 

6' = (0 2 ) to be independent from Ra and Pr, see Fig- 
ure 2J That figures shows that we have 9' ~ A for all 
Ra and Pr within our numerical precision. In contrast, 
Ref. @ predicted a dependence of the thermal fluctua- 
tions on both Ra and Pr, namely 0' / A ~ {PrRa)^ 1 / 8 for 
the regimes IV/ and IVJ which correspond to the bulk of 
turbulence analysed here. Our interpretation of Figure 
01 is that the bulk turbulence only has one temperature 
scale, namely A. For real RB turbulence it is the bound- 
ary layer dynamics which introduces further temperature 
scales, leading to the Ra and Pr number dependence 
of the temperature fluctuations observed in experiments 
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FIG. 3: 

a slope of 2.77 ± 0.03, slope 3 is shown for comparison, (b) 
ee I (kA 2 L -2 ) vs. RePr. We obtained a fitted slope 0.87±0. 04 
while slopes 1 and 1/2 are also shown for comparison. 



V. DYNAMICS OF THE FLOW 

In this section we provide an insight into the dynamics 
of the periodic Rayleigh-Benard flow. A bi-dimensional 
vertical snapshot of the flow is shown in Figure |S] Al- 
ready from this pictorial view the presence of an upward 
moving hot column and a downward moving cold column 
is clearly evident. 

Indeed these large scale structure can be related to the 
presence of "elevator modes" (or jets, forming in the flow) 
growing in time until finally breaking down due to some 
instability mechanisms. 

As proposed by Doering and collaborators in Q it 
is possible to predict the presence of these modes di- 
rectly starting from equations and (JHJ). Doering et 
dl. showed that, due to the periodic boundary condi- 
tions, this coupled system of equations admits a partic- 
ular solution 9 = 8oe KXt sin(fc ■ x), u^ — u$e KXt sin(fe • x), 
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FIG. 4: (a) Normalized temperature variance 9' / A vs. Ra 
at fixed Pr = 1. (b) 9' /A vs. Pr at fixed Ra = 1.4 ■ 10 7 . 



U2 = Ui = 0, which is independent from the vertical 
coordinate z (here k = [k x , k y )) and with: 



1 



(Pr + l)k 2 + 



(15) 



+ i,/(Pr + l) 2 fc 4 + 4Pr fe-fc 4 



From equation l|15|) one finds that the first unsta- 
ble mode appears for Ra > Ra c = (2tt) 4 ~ 1558.54, 
corresponding to the instability of the smallest possible 
wavenumber in the system, i.e. fc 2 = (2it / L) 2 n 2 with 
n=(l,0). 

The presence of accelerating modes with growth rate 
controlled by A can also be seen from Figure where 
we show Nu(t) on log-scale (notice the huge range over 
which Nu fluctuates). 

In Figure we show the PDF of Nu(t) which is 
strongly skewed towards large Nu values. This asym- 
metry reflects the periods of exponential growth (also 
visible in Figure As can be seen in Figure 00 for all 
Ra and Pr the system typically spends 54% of the time 
in growing modes. 

Also the relative fluctuations of Nu on the Ra and Pr 
numbers (see Figure EJ) seems to indicate no dependen- 
cies, at least in the range of parameters studied. 

Despite the presence of exact exploding solutions, our 
system clearly shows that in the turbulent regime these 
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FIG. 5: Snapshot of the flow, showing elevator modes and 
jets. Here 9 is shown in colors: red and yellow encode for 
positive values, with red greater in amplitude than yellow, 
green is for small values around zero, while blue stands for the 
negative values, the dark blue stands for the more negative 
values. Velocity in the same plane is shown with arrows. 



solutions become unstable due to some yet to be explored 
instability mechanism. The interplay between exploding 
modes and destabilization sets the value of the Nusselt 
number, i.e. the heat transfer through the cell. 

We stress that the study of the dynamics of the explo- 
sive solutions and of their successive collaps in a turbu- 
lent cell is crucial for the understanding the behaviour of 
"integral" quantities, like, for example, the heat transfer. 



VI. CONCLUSIONS 
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FIG. 6: (a) Time series Nu(t) for Ra = 1.4 • 10 7 (top) and 
Ra — 9.6T0 4 (bottom), in both case Pr — 1. (b) Logarithmic 
derivative of Nu(t) for Ra — 9.6T0 4 , here reproduced only for 
a small time section of the data in (a) . The series of horizontal 
lines represent the exponential rate of growing respectively 
(top to bottom) for the mode A(0, 1), A(0,2), and A(l,2). 



In conclusion, we confirmed that both the Ra- and 
the Pr-scaling of Nu and Re in homogeneous Rayleigh- 
Benard convection is consistent with the suggested scal- 
ing laws of the Grossmann-Lohse theory for the bulk- 
dominated regime (regime IV/ of |a Ilia), which is the 
so-called "ultimate regime of thermal convection". We 
also showed that the thermal and kinetic dissipations 
scale roughly as assumed in that theory. The temper- 
ature fluctuations do not show any Ra or Pr dependence 
for homogeneous Rayleigh-Benard convection. From the 
dynamics the heat transport and flow visualizations we 
identify "elevator modes" which are brought into the con- 
text of a recent analytical finding by Doering et al. In 
future work we plan to further elucidate the flow orga- 
nization and in particular the instability mechanisms of 
the elevator modes which set the Nusselt number in ho- 
mogeneous RB flow and therefore presumably also in the 



ultimate regime of thermal convection. 
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